COHERENT PIONS FROM NEUTRINO SCATTERING OFF NUCLEI 
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We describe a model for pion production off nucleons and coherent pions from nuclei induced 
by neutrinos in the 1 GeV energy regime. Besides the dominant A pole contribution, it takes into 
account the effect of background terms required by chiral symmetry. Moreover, the model uses 
a reduced nucleon-to-A resonance axial coupling, which leads to coherent pion production cross 
sections around a factor two smaller than most of the previous theoretical estimates. Nuclear effects 
like medium corrections on the A propagator and final pion distortion are included. 
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I. INTRODUCTION 

Neutrinos have been in the forefront of research in particle and nuclear physics for a long time. One of these fields 
is the study of pion production off nuclei induced by neutrinos. A proper understanding of this process is necessary in 
the analysis of the present generation of precision neutrino oscillation experiments. For instance, the 7r° produced in 
neutral currents (NC) is the most important z^-induced background to experiments like MiniBoone[l[ that are trying 
to measure —¥ v e oscillations in the neutrino energy range around 1 GeV. Also of importance is the background 
that appears from n + charged current (CC) production in — > v x disappearance searches like T2K0. Moreover the 
pion is strongly coupled to the A(1232) resonance, and neutrino scattering is presently the best way to access to the 
axial nucleon-A transition couplings. The most complete information in this regard comes from the bubble chamber 
data of ANL|3|, |4| and BNL 5, 6] where the target was cooled deuterium. However, in present oscillation experiments 
the target for neutrino interaction is finite nuclei, for instance 12 C (mineral oil in MiniBoone) or 16 (water target 
in T2K). This introduces sizable many-body effects that are difficult to disentangle from the genuine single nucleon 
response to the neutrino probe. 

In Sec. |TT]of these paper we describe a phenomenological model[7| for pion production induced by neutrino scattering 
off free nucleons. This model takes into account non-resonant background processes that are usually neglected. These 
non-resonant processes are determined by chiral symmetry and thus do not introduce free parameters. We then 
perform a fit of the axial A-A parameters and discuss a possible violation of the non-diagonal Goldbcrger-Treiman 
relation. In Sec. IIIII we extend the model to describe the CC coherent process in which the final nucleus is left in its 
ground state. We will further try to discuss how the coherent reaction can show some light on the values of the axial 
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iV-A couplings. We shall focus on the CC process, but the model can be easily extended to antineutrino reactions 
and NC processes. The interested reader could find further details in Refs. |7Hllj 

II. SINGLE NUCLEON PION PRODUCTION 

Here we review the model for the free nucleon reaction 

n(k)+ N^r(k')+N'+ tt+(£v) (1) 

as introduced in Ref. This model considers the dominant A pole mechanism in which the neutrino excites a 
A(1232) resonance that subsequently decays into Nir. In our model we have also included non-resonant background 
terms as required by chiral symmetry, see Fig. [TJ Some previous works |12Hl4j also considered background terms, 
though they were not consistent with the chiral counting. 

The vector part of the interaction N-A can be related to the electromagnetic current by imposing conservation of 
the vector current. For photon induced reactions extensive experimental data exist and in Ref. [15| they were employed 
to fit the vector current couplings. We shall use this fit in our work. Unfortunately the axial N-A current is not so 
well studied. The usual approach is to parameterize the interaction in terms of four form factors C% 4 56 (q 2 ). One 
can assume partial conservation of the axial curre nt ( PCAC) and obtain the relation Cq = C^M 2 /(m 2 — q 2 ), with 
M the nucleon mass. Furthermore, one can deduce |16j from dispersion relations the following conditions: C^(q 2 ) = 
and C± = — Thus we are left with only a free form factor, the dominant one C^(q 2 ). A different number of 

parameterizations have been proposed for this form factor, nevertheless the experimental data are quite limited and 
thus a simple dipole form 

should be enough. In order to keep the axial transition radius in the range 0.7-0.8 fm one expect the axial mass 
to have a value of around Ma ~ 0.85-1.0 GeV. Furthermore one can assume the well known Goldberger-Treiman 
relation (GTR) for ttNN coupling to be also valid for the nNA coupling and thus obtain 

C£(P) = J\f«— = 1-2, (3) 
V 3 m-x 

where f„ = 93MeV is the pion decay constant and /* = 2.2, the nNA coupling. Unfortunately, there are no 
constraints from Chiral Perturbation Theory [17j and the lattice QCD calculations fl8j are still inconclusive. 

Most of the approaches in the literature assume A dominance, that is, only include the first two diagrams in Fig. [T] 
We improve this situation by including non-resonant contributions 7] required by chiral symmetry. In addition to 
the A(1232) pole (AP) (first row) mechanism the model includes background terms required by chiral symmetry: 
nucleon (second row) pole terms contact and pion pole contribution (third row) and pion-in-flight term. We calculate 
them by using the SU(2) non-linear a model Lagrangian. The only parameters in the theory are the pion and nucleon 
masses and the pion decay constant. All other couplings are completely fixed by the theory, so no new parameter is 
introduced. 

We found (see Fig. [2]) that these background terms produced significant effects in all channels, namely an enhance- 
ment of about 10% in the cross section that resulted in a disagreement with the ANL data. As a result we had to 
readjust the strength of the dominant A pole contribution. The least known ingredients of the model are the axial 
nucleon-to-A transition form factors, of which gives the largest contribution. This strongly suggested a refit of 
that form factor to the experimental data, which we did by fitting the flux-averaged v^p — > n~pir + ANL ^-differential 
cross section for pion- nucleon invariant masses 32] W < 1.4 GeV. The obtained parameters were 

C^(0) = 0.87 ±0.08, Ma — 0.985 ± 0.082 GeV (4) 

with a x 2 /dof = 0.4 and a correlation coefficient r = —0.85, that amounts to a 30% reduction of the GTR prediction. 
Thus, our full model leads to an overall better description of the data for one-pion production reactions off the nucleon. 
This reduction of the C$(0) value is consistent with recent results in lattice QCD[l8| and quark models [l9|. 

Recently other fits have been proposed. For instance, in Ref. [2(| they keep the GTR but introduce a non-dipole 
form factor with additional parameters. As in the ANL data the relevant phase space is around q 2 =0.1 GeV 2 , they 
could keep the GTR at the cost of having a large dependence on q 2 for the form factor, that yields a large, somehow 
unphysical, axial transition radius of around 1.4 fm. Furthermore neither statistical errors nor correlation factors are 
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FIG. 1: Set of diagrams for the model for the W + N — > N'-k reaction. 




FIG. 2: Flux averaged q 2 — differential u^p — > [i~pn + cross section for the ANL (left) and BNL (right). Dashed lines stand for the 
contribution of the AP mechanism with the GTR assumption for C^. We also plot results with the full model of Fig.fl] assuming GTR 
(dashed-dotted) and with our best fit parameters, Eq. J4}. 



given in that reference. Another analysis plj raised new questions, namely the effect of deuterium wave function on 
the cross section and the flux uncertainties in the ANL and BNL data. The authors of this latter work took into 
account both effects, though we believe that their statistical analysis is not quite robust (see discussion in Ref. [llj). 
Furthermore they only took into account the dominant A contribution. The inclusion of deuterium wave function 
reduces the cross section about an 8%, so this somehow compensates the neglect of the non- resonant background, and 
they obtained a best fit of C^(0) = 1.19 ± 0.08 in agreement with the GTR assumption. 

Recently pd| we have improved our fit of Ref. 0, improving the lines suggested in Graczyk et al. In summary in 
this new fit 

1: all diagrams in Fig. Q] are included; 

2: the fitted data are the full ANL data set and the BNL total cross sections at the three lowest neutrino energies 
(we neglect higher energies where the effect of higher resonances beyond the A(1232) must be addressed; the 
BNL (^-differential cross sections were not taken into account as they are not normalized); 

3: deuterium wave function effects were introduced following the prescription of Ref. [22j |; 

4: the form factors C% and G\ were tentatively included in one fit, though the data were found to be quite 
insensitive to their values so we decided to stick to the Adler's assumption; and 

5: the uncertainties in ANL and BNL flux normalization are introduced as fully correlated systematic errors. 
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FIG. 3: Pion momentum (right) and angular (left) differential cross section. 



In this way we obtain a best fit of C5 (0) = 1.00 ± 0.11 and Ma = 0.93 ± 0.07 GeV with a goodness of fit value of 
X 2 d.o.f. = 0.42. Thus we observe a violations of the off diagonal Goldberger-Treiman relation at the level of 2a. 



III. THE COHERENT REACTION 

Here we describe our modelQ for the coherent reaction 

vi(k) + A z \ gs ( PA ) -> r(fc') + A z \ gs {p' A ) + 7r + (Av) (5) 

where the target nucleus Az is left in the ground state (gs). To calculate the amplitude of this process we sum over all 
individual nucleon wave functions, which are modelled by a Fermi gas in local density approximation. The individual 
nucleon amplitudes are modelled following the model of the previous section, using the fit of Eq. |4] for the form 
factor. On to p of that a number of many-body effects are introduced. In first place we take into account the in-medium 
modifications [23[ of the A(1232) properties. This implies a shift in the pole mass towards lower energies and most 
importantly a net broadening of the width (the opening of new decay channels in the nuclear medium compensates the 
Pauli blocking of the ttN decay channel). Also important is the distortion of the outgoing pion by strong interaction 
with the nucleus. Thus we consider the wave function of the pion to be the outgoing solution to the Klein-Gordon 



equation with a microscopic optical potential 24] whose imaginary part takes into account the inelastic interactions of 
the pion with the nucleus, that thus disappear from the coherent channel. We must emphasize here that solving the 
Klein-Gordon equation is the correct way of describing the distortion of the outgoing pion. Other approaches use either 
a Monte Carlo simulation [2^] or include an attenuation factor fitted to the pion nuclei scattering cross section (26|. 
The first procedure, though physically sound, can be a bit misleading as it includes in its cross section processes (like 
quasi-elastic scattered pions) that do not leave the nucleus in it ground state, thus are not coherent. This kind of 
models are used in the analysis of MiniBoonc experiment, thus making a bit messy the direct comparison between 
theoretical models and experimental results (see discussion in Ref. Q). The second approach is an oversimplification, 
as the pion-nucleus scattering is quite a different process from the neutrino induced pion production. As pion nucleus 
interaction is governed by the strong interaction, the incoming pion interacts strongly with the nuclear surface, thus 
the pion nucleus cross section is quite insensitive to the details of the nuclear core. On the other hand, neutrino 
scattering is a weak process, dominated by nuclear density, so pions are mostly produced in the deep, high density 
regions of the nucleus. The physics of pion interaction is thus quite different in pion scattering off nuclei and pion 
production by electroweak probes. 

In left panel of Fig. [3] we show the pion momentum distribution for CC coherent pion production, in the peak energy 
region of the T2K experiment. Including A in-medium self-energy (long-dashed line) reduces the PWIA results (short- 
dashed line). Further inclusion of pion distortion (full model, solid line) reduces the cross section, and the peak is 
shifted towards lower energies, reflecting the strong absorption and the higher probability of a quasi-elastic collisions 
of the pion in the A kinematical region. The total cross section reduction is around 60%. Similar nuclear effects 
were already studied in Refs. [27l |28j |. However, the authors of these references neglected the nucleon momenta in 
the Dirac spinors. The effect of this approximation (nucleons at rest, dotted line) results in a ~ 15% decrease of the 
total cross section. In the right panel of Fig. [3] we show the pion angular distribution with respect to the incoming 
neutrino direction. The reaction is very forward peaked, as expected due to the nucleus form factor. The angular 
distribution profile keeps its forward peaked behavior after introduction of nuclear medium effects. Furthermore we 
corrected some numerical errors in the mentioned papers. However one must be aware that this model does not take 
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Reaction Experiment a(10 40 cm 2 ) a(10 40 cm 2 ) Exp 



cc 


^+ 12 c 


K2K 




NC 


^+ 12 c 


MiniBoone 


3.33 


CC 


^+ 12 c 


MiniBoone 


4.46 


cc 


^+ 16 o 


T2K 


4.19 


cc 




T2K 


3.54 



TABLE I: Total cross sections for the coherent process. We neglect the highest 10% of the energy spectrum. 



into account the non-localities [29j in the A propagation. We believe this effect is partially taken into account in an 
effective fashion by our treatment of the A in nuclear medium; nevertheless further studies would be interesting. 

In Table Mil we compare our model with present results of K2K experiment [3(| and show our predicitions for 
MiniBoone and T2K. Our prediction, subject to sizable uncertainties, lies well below the K2K upper bound, due to 
the use of a low value for C^(0), while our prediction for the NC MiniBoone cross section is notably smaller than 
that given in the PhD thesis of J. L. Raaf[3l|. However, we believe (see discussion in Refs. d,[l3 that the MiniBoone 
analysis might importantly overestimate this cross section, not only because some of the ir°s which undergo FSI 
collisions are accounted for instead of being removed, but also because a possible mis-match between the absolute 
normalisation of the background and coherent yields. Note that the K2K and MiniBoone results seems somehow 
incompatible with the approximate relation ctqc ~ 2ctnc , which would be expected from A dominance and neglecting 
finite muon mass effects. 
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